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ABSTRACT 

We present a study of outflow, infall, and rotation in a ~10 5 L & star-forming 
region, IRAS 18360-0537, with Submillimeter Array (SMA) and IRAM 30m ob- 
servations. The 1.3 mm continuum map shows a 0.5 pc dust ridge, of which 
the central compact part has a mass of ~80 M and harbors two condensations, 
MM1 and MM2. The CO (2-1) and SiO (5-4) maps reveal a biconical outflow 
centered at MM1, which is a hot molecular core (HMC) with a gas temperature 
of 320 ± 50 K and a mass of ~13 M & . The outflow has a gas mass of 54 M & and 
a dynamical timescale of 8 x 10 3 yr. The kinematics of the HMC is probed by 
high-excitation CH 3 OH and CH 3 CN lines, which are detected at sub-arcsecond 
resolution and unveil a velocity gradient perpendicular to the outflow axis, sug- 
gesting a disk-like rotation of the HMC. An infalling envelope around the HMC 
is evidenced by CN lines exhibiting a profound inverse P-Cygni profile, and the 
estimated mass infall rate, 1.5 x 10~ 3 M & yr" 1 , is well comparable to that inferred 
from the mass outflow rate. A more detailed investigation of the kinematics of 
the dense gas around the HMC is obtained from the 13 CO and C 18 (2-1) lines; 
the position-velocity diagrams of the two lines are consistent with the model of a 
free-falling and Keplerian-like rotating envelope. The observations suggest that 
the protostar of a current mass ~10 M embedded within MM1 will develop into 
an O star via disk accretion and envelope infall. 

Subject headings: stars: formation — stars: early-type — ISM: kinematics and 
dynamics — ISM: jets and outflows 
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Introduction 



The formation process of high-mass stars (> 8 M Q ) is poorly understood ( IZinnecker fc Yorke 20071 ) 
One key question in the field is whether high-mas s stars form as a sca led-up version of 



low-mass star formation (e.g., iKeto fc Zhang 201CH ; I Johnston et al. 201 if ). In a standard 



paradigm of low-mass star fo rmation, a protos tar accretes from a circumstellar disk embedded 
within an infalling envelope ( jShu et al. 1987T ) . Since accretion disks and collimated outflows 
are physically connected, the detection of well collimated outflows in an increasing number 
of high-mass star-forming regions supports an accretion-based scenario for massive star for 



mation (e.g.. IBeuther et al. 2002c IZhang et al. 20 07a; 



Qiu fc Zhang 20091 ; iQiu et al. 2011aT IZapata et al. 2011 



iang et al. 2007bl : biu et al. 2007 



In the meantime, great efforts 
have been made in searching for direct evidence for an infalling envelope and a rotating disk in 
association with massive star formation. Nevertheless, only a very few massive star-forming 
cores have been found to show clear signatures of both infall and rotation. Beltran (2011) en- 
lists a sample of six such sources from the literature, G10.62-0.38, G19.61-0.23, G24.78+0.08 
Al, G31.41+0.31, W51 e2, and W51 North. But some more cases exist (e.g., NGC 7538 



IRS 1, IQiu et al. 201 lbl ; IBeuther et al. 20121 ). In these studies, a rotation is inferred from 
a velocity gradient seen in high-density tracing molecular lines and an infall is probed by 
molecular lines with redshifted absorption against a bright continuum source. Even this 
small sample should be taken with caution. For example, the orientations of the velocity 
gradients may change by tens of degre e s, depending on the lines being analyzed (e.g., W51 e2: 
Zhang fc Ho 19971 Izhang et al. 19981 : IKeto fc Klaassen 20081 : iKlaassen et al. 20091 ). Mean- 
while, if a velocity gradient is arising from a rotating disk, it is expected to be perpendicular 
to a bipolar outflow, but for most sources in the sample, such an as sociation between a p uta- 
tive disk and an outflow is yet to be confirmed (e.g., G31.41+0.31. ICesaroni et al. 201ll and 
references therein). Therefore, identifying and characterizing clear-cut examples showing 
self-consistent outflow, infall, and rotation motions remains of great interest to our under- 
standing of massive star formation. Here we present such a study toward IRAS 18360-0537, 
using Submillimeter Arra;y0 (SMA) 0."5-l."5 resolution observations in the 1.3 mm waveband. 

IRA S 18360-0537 has a k inematic distance of 6.3 kpc and a far-IR luminosity of 1.2 x 
10 5 L Q (IMolinari et al. 19961 ). In earli er single-dish surveys, the region wa s found to be 
associ ated with water ma ser emissio n (IPalla et al. 199ll: [Brand et al. 1994 ) and thermal 
NH 3 (IMolinari et al. 19961 ) and CS (IBronfman et al. 19961 ) emission. There is so far no 
report on high-angular-resolution observations of the region in (sub)millimeter or centimeter 



6 The SMA is joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia 
Sinica. 
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wavelengths. 



2. Observations 

The SMA observations were undertaken in two frequency setups, each with two array 
configurations. The observing campaign was designed to cover CO, 13 CO, C 18 (2-1), and 
several CH3OH and CH3CN lines in one frequency setup and to cover CN (iV=2-l) and 
SiO (5-4) lines in the other setup. The CO and SiO lines are typical outflow tracers, 
while the CN lines may reveal pote ntial infall signatures in high-mass star-forming cores 



(jZapata et al. 20081 ; IWu et al. 20091 ). The kinematics of the innermost dense gas can be 
inferred from high-excitation CH3OH and CH3CN lines. The details of the observations are 
presented in Table [TJ 

The data were calibrated using the IDL MIR package, and imaged with the MIRIAD 
software. A continuum data set was constructed from line free channels using the MIRIAD 
task UVLIN. We performed self-calibration on the continuum, and applied the solutions to 
both continuum and line data sets. For molecular line maps within each frequency setup, 
the calibrated visibilities from the two array configurations were jointly imaged. The angular 
resolution (the synthesized beam, 9 syn , and its position angle, PA) and rms noise level, a, of a 
map vary with the visibility weighting algorithm. Most maps were made with a compromise 
between resolution and sensitivity (e.g., varying the ROBUST parameter), resulting in 9 syn ~ 
2."7 x l."4 at PA - 62°, a ~ 50 mJybeam -1 per 1.5 kms -1 in CO lines, 6 syn - 2."0 x l."4 
at PA ~ 64°, a ~ 30 mJybeam -1 per 1.5 kms -1 in SiO (5-4), and 9 syn ~ l."4 x l."3 at 
PA ~ 57°, o ~ 35 mJybeam -1 per 0.6 kms -1 in CN lines. To achieve the highest possible 
angular resolution, which is desirable in searching for disk-like structures, the CH 3 OH and 
CH3CN lines were made with a uniform weighting, leading to 8 syn ~ 0."8 x 0."6 at PA ~ 20°, 
a ~ 65 mJybeam -1 per 0.6 kms -1 . For the 1.3 mm continuum emission, we made a map 
combining visibilities from all four array configurations (Figure [T]( a)), with 8 syQ ~ l."6 x l."2 
at PA ~ 68°, o ~ 1.5 mJybeam -1 , and a map combining the EXT and VEX data (Figure 
E»), with 9 syn ~ 0."6 x 0."5 at PA ~ 18°, and a ~ 2.7 mJybeam -1 . 



The IRAMLj 30m telescope was used to fill the zero-spacing of the SMA maps in CO 
and 13 CO (2-1). The observations were conducted on 2007 November 1. A 2' x 2' region 
centered at (R.A., Dec.)j 2 ooo=(18 h 38 m 40. s 3, — 5°35'6."0) was mapped in the on-the-fly mode. 
Scans were made in both Right Ascension and Declination directions in an effort to prevent 
systematic scanning effects from appearing in the data. The receiver was tuned to 230.538 



7 IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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GHz and 220.399 GHz for the CO and 13 CO (2-1) lines, respectively, and the spectra have 
a 0.4 kms -1 resolution. The spectra were processed using CLASS in the GILDAS software 
package. The CO spectra have an rms noise level of 1.1 K and the 13 CO spectra have an 
rms of 0.5 K in T m b- The combination of the SMA and single-dish data was performed in 
MIRIAD following a procedure outlined in Zhang et al. (1995). 



3. Results 

3.1. Dust Condensations 

Figure EJa) shows the 1.3 mm continuum map made from observations conducted with 
the four array configurations. The emission reveals a northwest-southeast (NW-SE) ridge 
with a bright and compact part at the center. In Figure [T^b), relatively extended emission 
is filtered out in the map made with the EXT and VEX data, and the central part is 
resolved into two condensations, MM1 and MM2. This work focuses on the brightest dust 
condensation, MM1. It shows very faint free-free emission (~0.5 mJy at 1-3 cm) in our Very 
Large Array observations (Qiu et al. in prep.). The observed 1.3 mm continuum emission is 
completely dominated by dust emission. 

MM1 shows rich emission lines in complex organic molecules, e.g., CH3OH, CH3CN, 
C2H3CN, and C 2 H 5 CN. The CH 3 CN (12-11) fT-ladder is a good thermometer of dense gas 
in high-mass star-forming regions. Following Qiu & Zhang (2009) and Qiu et al. (2011b), we 
simultaneously fit the K = 2-8 components^], taking into account the optical depth effect. 
The systemic velocity, 103.5 kms -1 , is determined from Gaussian fittings to the K — 7, 8 
components, which are of the highest excitation and not heavily blended with other lines. 
This determination is concurred by the C 2 H 5 CN (252 2 ,4-24 2 2,3) and CHj 3 CN (12 2 -11 2 ) lines 
(see Figure H]), to which Gaussian fittings give velocities of 102.9 and 103.5 kms -1 , respec- 
tively. In Figure HJ the best-fit model reasonably matches the observation, and yields a gas 
temperature of 320 ±50 K, where the uncertainty accounts for the la noise level. The fitting 
assumes that all the K components are tracing the same gas under local thermodynamic 
equilibrium conditions. Provided potential density and temperature gradients for an inter- 
nally heated protostellar core, the best-fit model may represent a characteristic estimate. 

Assuming thermal equilibrium between the dust and dense gas, the derived gas tem- 
perature in MM1 approximates the dust temperature. To constrain the dust opacity index, 
0, we make use of the 880 fim continuum data obtained from recent SMA observations. 



The relatively low excitation K = 0, 1 components are resolved out in the VEX data. 
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By comparing the 1.3 mm and 880 /zm continuum maps constructed with the same (u,v) 
range, we deduce (3 = 0.82 toward MM1 and (3 = 0.50 to ward MM2, corres ponding to dust 
opacities of 2.5 and 4.3 cm 2 g _1 , respectively, at 1.3 mm ( jHildebrand 19831 ). The gas mass 
of MM1 is then estimated from its 1.3 mm dust continuum flux, 0.76 Jy, to be 13 M Q , by 
adopting a canonical gas-to-dust mass ratio of 100. The mass of MM2 is less clear given 
its unknown gas o r dust temperature. A dopting a temperature of 30 K deduced from prior 
NH 3 observations (IMolinari et al. 19961 ). one obtains a mass of 25 M from the continuum 
flux of 0.2 Jy. 



3.2. Bipolar Outflow 

Figure |3] shows the velocity integrated CO (2-1) and SiO (5-4) emission observed with 
the SMA. Both CO and SiO emission seems to be dominated by a northeast-southwest 
(NE-SW) bipolar outflow centered at MM1. The outflow has a biconical shape with an 
opening angle of ~65°. Such an outflow is often thought to form from ambient gas b eing 



swept up or entrained by a wide-angle wind (e.g., IShepherd et al. 19981 ; lOiu et al. 20091 ). A 
blueshifted clump to the south of MM1 and a redshifted clump to the north, which are more 
prominent in the CO map, could be part of the biconical outflow, if the wide-angle wind is 
slightly inclined to the plane of sky and impinges into an environment with an inhomogeneous 
density structure. However, the possibility that the two clumps form another outflow cannot 
be ruled out. 

Figure H] shows the velocity channel maps of the CO emission observed with the SMA; 
at low velocities (96-112.5 kins" 1 ), the maps made with the combined SMA and IRAM 
30m (SMA+30m) data are also presented. It is clear that the combined data effectively 
recover extended emission around the cloud velocity. Emission arising from outflow struc- 
tures appears predominant in channels of <100.5 kms -1 and >106.5 kms -1 . In particular, 
the cone-shaped structure of the SW lobe is most noticeable in the 106.5 kms -1 channel 
map. The emission at 111.0 and 112.5 kms -1 is mostly filtered out in the SMA maps and 
appears to be dominated by diffuse gas in the combined maps. This is most likely due to the 
self-absorption effect caused by a foreground cloud. Figure |5] shows the SiO (5-4) channel 
maps. SiO is believed to be produ ced in the gas phas e through shock-driven sputtering of 



Si-bearing species on dust grains (jSchilke et al. 19971 ) . Therefore the SiO emission is not 



contaminated by emission from quiescent ambient gas. As seen in the integrated map, the 
SiO outflow is in general consistent with the outflow seen in the CO emission. 

Figure |6] shows the SMA 13 CO (2-1) channel maps as well as the SMA+30m maps at 
low velocities. The 13 CO emission arising from the NE-SW outflow is less extended than 
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that in CO (2-1), but is still appreciable in the channels of 94.5-108.0 kms -1 (except the 
103.5 kms -1 channel). In the close vicinity of MM1, the 13 CO emission peaks to the NW of 
MM1 at blueshifted velocities and to the SE at redshifted velocities; the offset between the 
13 CO peak and dust continuum peak decreases with the velocity so that the 13 CO emission 
approximately coincides with MM1 at the highest velocities. On the other hand, the CO 
(2-1) emission peak is clearly offset from MM1 even at the highest velocities, and shifts from 
the NE to SW. Such a difference is highlighted in Figure [7J While the relatively extended 
emission in 13 CO traces the outflow, the compact and bright emission closely around MM1 
is predominantly arising from a dense envelope. 



3.3. Velocity Gradients in High-excitation CH 3 OH and CH 3 CN Lines 

The 13 CO (2-1) observations provide hints for the rotation of dense gas around MM1 
(Figure [7J. However, the 13 CO (2-1) line is affected by the outflow and, with a critical 
density, n cr , of order 10 4 cm -3 and an upper level energy, E up , around 16 K above the ground, 
cannot effectively probe the hot and dense gas within MM1. To search for a potential disk- 
like structure, which is supposed to lie in the innermost part of a molecular core, we inspected 
high-excitation lines from less abundant molecular species. A NW-SE velocity gradient across 
MM1 is consistently seen. Figure El shows the zeroth moment (velocity integrated emission) 
and first moment (intensity weighted velocity) maps, and position-velocity (P-V) diagrams 
of the CH 3 OH (154,11-163,13)^, CH 3 CN (12 7 -11 7 ), and vibrational CH 3 CN (121-110 lines. 
These lines are chosen for the following reasons. First, they are highly excited {E up > 350 
K) and, as evidenced by the zeroth moment maps, probe the innermost dense gas in MM1. 
Second, they are not significantly blended with other lines; CH3CN (I27-II7) is slightly 
blended with CH^CN (12 5 -11 5 ), but the two lines are ~9.5 kms -1 apart and the former is 
three times brighter than the latter, so it is feasible to ignore the contribution from CH^CN 
(I25-II5) in interpreting CH3CN (I27-II7). Finally, these lines are detected with sufficient 
signal-to-noise ratios in the VEX observations so that the velocity gradient is discernable at 
a 0."4 resolution; some other lines with a higher E np , e.g., CH 3 CN (12 8 -11 8 ), or with a lower 
optical depth, e.g., CH^CN (I22-II2), are subject to limited signal-to-noise ratios and barely 
show the NW-SE velocity gradient. From the first moment maps of CH3OH (154,11-163,13)^ 
and CH3CN (12 7 -11 7 ), the velocity gradient is approximately along a position angle (PA) of 
140°, which is adopted for the cut of the P-V diagrams. The vibrational CH 3 CN (12]— 110 
line has the highest excitation and lowest signal-to-noise ratio among the three lines shown 
here. The direction of the velocity gradient in its first-moment map is less clear, but an 
overall NW-SE gradient is still discernable and the P-V pattern is consistent with those 
of the other two lines. Compared to the orientation of the biconical NE-SW outflow, the 
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CH 3 OH and CH 3 CN velocity gradient is approximately perpendicular to the outflow axis, 
and is apparently arising from a disk-like rotation. However, neither the dust continuum nor 
the molecular line maps of the MM1 condensation shows a flattened morphology. This is 
very likely due to insufficient angular resolutions (> 0."5) of the observations. 



3.4. Infall Signatures Seen in CN (JV=2-1) 

Several CN (N=2-l) hyperfine lines are detected toward MM1. The brightest three, 
i.e., the J=5/2-3/2, F=5/2-3/2 (z/ =226.87417 GHz), J=5/2-3/2, F=7/2-5/2 (226.87474 
GHz), and J=5/2-3/2, F=3/2-l/2 (226.87590 GHz) lines, are heavily blended and form 
a triplet. Besides the triplet, the lines of J=3/2-l/2, F=5/2-3/2 (226.65956 GHz) and 
J=3/2-l/2, F=3/2-l/2 (226.67931 GHz) are detected as well; these two lines do not suffer 
from line blending and the brighter one, CN (N=2-l, J=3/2— 1/2, F=5/2-3/2), is shown 
in Figure [9J The CN (JV=2— 1) transitions have E np ~ 16 K and n cr of order 10 7 cm" 3 , 
probing relatively cold gas in a dense envelope around MM1. All the detected CN lines show 
emission at blueshifted velocities and absorption at redshifted velocities, i.e., an inverse P- 
Cygni profile. Such a profile is often interpreted to indicate an infall of matter, whereas 
the absorption mostly comes from the gas on the ne ar side of a centra l source and the 



emission is attrib uted to the gas in the far side (e.g., IZapata et al. 20081 ; IWu et al. 2009 



Qiu et al. 201 lbl ). Therefore, the CN lines probe an infalling envelope around MM1. In 
Figure [9j the peaking velocities of the emission and absorption are slightly asymmetric with 
respective to the systemic velocity. It could be in part due to the uncertainty in determining 
the systemic velocity, which is estimated to be 103.5 kms " 1 in our data but w as suggested 



to be 102.3 kms 1 based on single-dish NH 3 observations (jMolinari et al. 19961 ). 



4. Discussion 

4.1. An accreting star of ~10 M Q 

MM1 has a gas temperature of 320 K and molecular chemistry characteristic of a hot 
molecular core (HMC). At least one high- mass protostar has formed in MM1. On the other 
hand, only very faint and spatially unresolved radio continuum emission is detected. The 
emission can be arising from an ionized jet or wind, or from a hypercompact H n region. In 
the latter case, the central ionizin g star has an equivale nt spectral type of B0. 5, corresponding 



to a stellar mass, M*, of 12 M (ISchaller et al. 19921 ). The mass of the outflowing gas can 



also shed light on the stellar mass of the protostar. Assuming that the molecular outflow 
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is momentum driven by an underlying wind, the accumulated mass ejected to the wind is 
M out i; out /t> w , where M out is the outflow mass, t> out the outflow velocity, and v w the wind 
speed. The mass loss rate to the wind is a fraction, /, of the mass infall rate, thus the mas s 



of the forming star can be expressed as M out v out /v w ■ (1 — /)// (see also lLada &: Fich 19961 ). 
The outflow mass is calculated from the CO line wings (<100.5 kms -1 and >106.5 kms -1 ), 
assuming an excitation temperature of 30 K and a CO abundance of 10 -4 . To correct for 
the optical depth effec t, we compare the CO and 13 CO (2-1) fluxes and adopt a C-to- 13 C 



abundance ratio of 37 (IWilson fc Rood 19941 ). A more detailed description of the calculation 
procedure can be found in Qiu et al. (2009). We derive an outflow mass of 54 M Q with the 
combined SMA and IRAM 30m data. The maximum velocity of the outflow observed in CO 



(2-1) is a bout 36 kms 1 . Adopting a wind speed of 500 kms 1 (e.g.. iMartf et al. 19981 ) and 



/ = 1/3 (ITomisaka 19981 ; IShu et al. 20001 ). the accumulated stellar mass is about 8 M , in 



agreement with the above estimate. 

The stellar luminosity of a ~10 M star at the Zero-Age-Main-Sequence (ZAMS) 
is ~10 4 Lq, only a small fraction of the total luminosity observed in IRAS 18360-0537 
(~10 5 Lq). The majority of the observed luminosity presumably comes from accretion 
shocks. The accretion luminosity can be estimated from GM^M / i?*, where M is the mass 
flux of the accretion onto the star and R* the radius of the star. Again, under the assumption 
of momentum conservation between the molecular outflow and the underlying wind, M can 
be inferred from the mass outflow rate. Without correcting for an unknown inclination angle, 
the dynamical timescale of the outflow, i.e., a half of the outflow extension (~0.3 pc) divided 
by the terminal velocity, is about 8 x 10 3 yr. The mass outflow rate is then 7 x 10 -3 M & yr -1 , 
leading to a mass infall rate of 1.5 x 10~ 3 M & yr -1 a nd M ~ 1.0 x 10 ~ 3 Mp, yr" 1 . If we adopt 



a radius of 4i? for a 10 M star at the ZAMS (IPanagia 19731 ; ISchaller et al. 19921 ). the 
accretion luminosity amounts to 8 x 10 4 L , which does appear to dominate the total lumi- 
nosity. 

Alternatively, according to the model of Hosokawa et al. (2010), massive protostars with 
disk accretion rates of 10 -3 M & yr -1 attain a total luminosity of ~10 5 L & as M* increases to 
~15 M Q , and reach the ZAMS for M* ~ 30 M . In this model, the stellar radius swells to 
tens of R & as M* reaches ~ 10 M , so the accretion luminosity is only of order 10 4 L , while 
the stellar radiation, which comes from the release of the gravitational energy, is dominating 
the total luminosity. The large radius of the protostar also reduces the effective temperature 
and the UV luminosity, thus the faint radio continuum emission is expected to be arising 
from an ionized jet or wind. 



Nevertheless, the protostar embedded in MM1 appears to have a current mass of ~10 
M & and is in an active accretion phase. 
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4.2. Rotation and Infall 

Is the accretion onto the central high-mass protostar mediated by a rotating disk? The 
observed CO and SiO outflow supports a positive answer since a biconi cal outflow is exp ected 



to originate from an accretion disk in low-mass star formation (see iQiu et al. 20091 . for a 
detailed comparison between wide-angle outflows in low-mass and high- mass star formation). 
The high-excitation CH3OH and CH3CN lines observed at sub-arcsecond resolution provide 
direct evidence for the rotation of the HMC There seem to be hints of a Keplerian-like 
motion in the P-V diagrams of the high-excitation lines. For example, in the upper right 
quadrant of Figure [S][e) , the maximum velocity increases with the decreasing offset; in the 
lower left quadrant of Figure E^d), a secondary peak around 98.5 kms -1 could be attributed 
to a Keplerian-like motion as well. 

A more detailed investigation of the rotation curve of the HMC awaits future observa- 
tions with greatly improved resolution and sensitivity. Here we look into the lines from more 
abundant molecular species, i.e., C 18 and 13 CO, which may help t o deliver a more com - 



plete view of the velocity field for the dense gas on a larger scale (e.g.. ICesaroni et al. 201 ih . 
13 CO (2-1) can be contaminated by the outflow, but from Figure [7J the emission in the close 
vicinity of MM1 is dominated by the dense envelope. In Figure ITUT a) we plot the 13 CO and 
C 18 P-V diagrams, cut along the PA of 140°. We note that the 13 CO and C 18 emission 
arises from a region encompassing both MM1 and MM2. However, from the P-V diagrams 
the velocity field cannot be ascribed to an orbiting motion of a binary. In particular, the 
high-excitation lines shown in Figure |8] are not detected toward MM2; these lines show a 
clear velocity gradient across MM1. In Figure [101 the CH 3 CN, C 18 0, and 13 CO emission 
seems to be arising from different layers of a coherent structure centered at MM1, whereas 
MM2 is likely a fragment of this structure. Despite the missing flux closely around the cloud 
velocity, the 13 CO P-V diagram shows a concave outer edge and is slightly skewed in a sense 
consistent with the velocity gradient seen in CH3OH and CH3CN. A similar pattern is dis- 
cernible in the C 18 P-V diagram, though the emission has a lower velocity dispersion and a 
smaller spatial extent compared to th e 13 CO emission. Such a P-V pattern ma y characterize 



a combination of infall and rotation ( ICesaroni et al. 201 ll ; iTobin et al. 20121 ). To corrobo- 



rate this interpretation, we overlay in Figure [TUl a free-fall and Keplerian-like rotation model 
which is described by Cesaroni et al. (2011). Since the C 18 and 13 CO emission mostly 
traces outer layers around the HMC, we set a central dynamical mass of 25 M , which is 
approximately the combined mass of the HMC and the embedded high-mass protostar, and 
an inner truncation radius of 0."25 (0.008 pc), which is approximately the deconvolved radius 
of the CH3OH and CH 3 CN emission; the outer radius is set to be 3."5 (0.1 pc), roughly the 
radius of the C 18 and 13 CO emission. The constructed model (solid curve) reasonably fits 
the C 18 and 13 CO P-V diagrams except the high velocity tail close to the zero offset seen 
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in 13 CO. That high velocity tail can be partly ascribed to the inner rotating and infalling 
gas, as illustrated by a scaled-down model (dashed curve in Figure [10]) with a central mass 
of 10 M and outer and inner radii of 0."25 and 0."032 (200 AU), respectively, but it can also 
be arising from the outflowing gas. Overall the bulk 13 CO and C 18 emission delineates an 
infalling and rotating envelope around the HMC 

To assess the envelope mass, we performed a two-component Gaussian fitting to the 
central compact part of the dust ridge with the MIRIAD task IMFIT, and derived a flux 0.96 
Jy for the component encompassing MM1. Subtracting the HMC contribution (0.76 Jy, see 
Section I37L]) . the dust continuum flux from the cool envelope amounts t o 0.2 Jy, resulting in a 



gas mass of 42 M for (3 = 0.82 and an adopted temperature of 30 K (IMolinari et al. 19961 ). 
From the fitting, the MM2 component has a flux of 0.2 Jy, consistent with the measurement 
based on Figure [T^b) (Section 13. ip . Hence, the total mass of the central compact part of 
the dust ridge reaches 80 M , whereas the mass of the envelope traced by the 13 CO and 
C 18 emission amounts to 67 M . Such a massive envelope is presumably gravitationally 
unstable. Indeed, MM2 is likely the consequence of the fragmentation of the envelope. 
However, the HMC itself has a gas mass (13 M ) comparable to the mass of the accreting 
protostar (10 M ), thus the latter can have an appreciable stabilizing effect on the former. 
Cesaroni et al. (2007) suggest that a rotational structure of a mass comparable to that of 
the central protostar, which is in turn less than 20 M , resembles an accretion disk rather 
than a transient toroid. In this sense, we may expect that as the central protostar continues 
to accrete and increase its mass, the rotating HMC serves as an ac cretion disk, and finally 



dissipates by enhanced UV radiation (e.g., by photoevaporation, iHollenbach et al. 19941 ) 
rather than by gravitational fragmentation. 

The infall motion of the envelope around the HMC is not only appreciable in the C 18 
and 13 CO P-V diagrams, but also more unambiguously evidenced by the CN hyperfine lines 
with a profound inverse P-Cygni profile. The mass infall rate can be estimated from 4irr 2 pVi n , 
where the infall velocity V\ n is about 1.5 kms -1 according to the CN (A=2-l) profiles (e.g., 
Figure [9]), r the radius of the envelope, and p the mass density. Considering that the CN 
lines were imaged with a synthesized beam (l."4 x l."3) similar to that of the dust continuum 
map in Figure Ufa) (l."6 x l."2), we estimate the radius and density based on the above two- 
component Gaussian fitting. For the component encompassing MM1, r ~ 0."9 (5700 AU) 
from the fitting, and the averaged number density is 8 x 10 6 cm~ 3 for a gas mass of 42 M . 
The mass infall rate is then found to be about 1.5 x 10~ 3 M yr _1 . From the model shown 
in Figure [TU| the free-fall velocity reaches 3 kms -1 at r ~ 0."9, a factor of two higher that 
inferred from the CN profiles. Such a discrepancy can be ascribed to the deviation of the 
infall motion from the modeled free-fall, the uncertainty of the systemic velocity (see Section 
13. 4p . or an underestimate of the radius from the two-component fitting. Nevertheless, since 
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the (free-fall plus Keplerian-like rotation) model parameters are not stringently constrained 
and the CN profiles provide an averaged estimate, the discrepancy is not significant. On 
the other hand, as discussed above, the mass outflow rate implies a mass infall rate of 
1.5 x 10~ 3 M yr _1 . Though both are subject to large uncertainties, the two independent 
estimates of the mass infall rate appear to agree with each other, suggesting that the envelope 
infall is feeding the accretion onto one protostar, i.e., the ~10 M© star embedded in MM1. 
Provided its high accretion rate and the massive gas reservoir available from the envelope, it 
is reasonable to expect that the star would significantly increase its mass and develop into 
an O star. 



5. Summary 

We have presented high-angular-resolution dust continuum and molecular line observa- 
tions toward IRAS 18360-0537, a ~10 5 L Q star-forming region. The region shows a 0.5 pc 
dust ridge, harboring two massive condensations, MM1 and MM2, in the central part. MM1, 
a HMC with a gas temperature of 320 K and a mass of 13 M , is embedded within an infalling 
envelope, and appears to be rotating about the axis of a bipolar outflow. 

The outflow is seen in CO (2-1) and SiO (5-4) and shows a biconical shape centered at 
MM1. The gas mass of the outflow is estimated to be 54 M , and the dynamical timescale 
is 8 x 10 3 yr. The rotation of the HMC is inferred from sub-arcsecond observations of 
high-excitation CH3OH and CH3CN lines, which show a velocity gradient approximately 
perpendicular to the outflow axis. On the other hand, the infall motion of the envelope 
around the HMC is probed by CN (iV=2-l) hyperfine lines, which exhibit a profound inverse 
P-Cygni profile. The estimated mass infall rate is 1.5 x 10~ 3 M Q yr _1 , comparable to that 
inferred from the mass outflow rate. Further information on the kinematics of the gas 
around the HMC is obtained from the C 18 and 13 CO (2-1) lines, whose P-V diagrams can 
be described with a model of a free-falling and Keplerian-like rotating envelope. With these 
observations, we speculate that a protostar in MM1, which we suggest to have a current 
mass of ~10 M Q , will grow into an O star via disk accretion and envelope infall, i.e., through 
a process similar to the standard paradigm for low-mass star formation. 
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Right Ascension (J2000) 

Fig. 1. — (a) 1.3 mm continuum map constructed by combining data from all the obser- 
vations. Solid and dotted contours show positive and negative emission, respectively, with 
the contour levels of 6 x (1, 2, 3, 23) L5 mJybeaur 1 . (b) Same as (a), but for the EXT 
and VEX data, and contour levels of 9 x (1,2,3, ...,21) 13 mJybeam -1 . In each panel (and 
hereafter), two filled squares denote the peak positions of the two condensations (MM1 and 
MM2), and a filled ellipse in the lower left shows the synthesized beam. 
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220.65 220.60 220.55 220.50 220.45 220.40 

Frequency (GHz) 

Fig. 2.— The spectra toward MM1, covering the K=2-8 ladder of CH 3 CN (12-11), with 
the observations shown in black histograms and the best fit LTE model in dashed lines. The 
CH 3 CN (12-11) K=2-8 components are labeled with solid bars, and the CH^CN (12-11) 
K=0-6 components with dotted bars. Some other lines covered in the frequency range are 
also denoted. 
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Right Ascension (J2000) Right Ascension (J2000) 

Fig. 3. — (a) CO (2-1) emission observed with the SMA, and integrated from 69 to 99 
kms -1 for the blueshifted lobe and from 106.5 to 139.5 kms -1 for the redshifted lobe, shown 
in blue and red contours, respectively. The starting and spacing contour levels are 10% 
of the peak emission, which is 56.4 Jybeam -1 kms -1 for the blueshifted lobe and 71.3 
Jybeam -1 kms -1 for the redshifted lobe. The gray curves outline a biconical outflow, and 
the two arrows indicate the orientation of the outflow axis. A dotted line marks the direction 
of a velocity gradient seen in MM1 (see Figure (b) The same as (a), but for SiO (5-4). 
The blueshifted lobe is integrated from 88.5 to 102 kms -1 , and the redshifted lobe from 105 
to 118.5 kms -1 . The starting and spacing contour levels are 10% of the peak emission, which 
is 6.7 Jybeam -1 kms -1 for the blueshifted lobe and 7.3 Jybeam -1 kms -1 for the redshifted 
lobe. 
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Fig. 4. — SMA CO (2-1) velocity channel maps. For channels of 96-112.5 kms -1 , the maps 
made with the combined SMA and IRAM 30m data are shown on the bottom. Solid and 
dotted contours show the positive and negative emission, respectively, with the same absolute 
levels starting from and increasing in steps of 0.6 Jybeam -1 . 
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Fig. 5. — SMA SiO (5-4) velocity channel maps. Solid and dotted contours show the 
positive and negative emission, respectively, with the same absolute levels starting from and 
increasing in steps of 0.12 Jybeam -1 . 



- 19 - 



'85:5 : 

eg, 


'a™' 


88:5 




'91.5 : 


CD- 
CO 

... 1 .... 1 .... 


:&5 : 


ko' ' J 




'Vos'b^"" 


:1 00.5' '1^^^ 

: «^ 


102.feS 


[\m ! 




_ ^_ 




'V69.5 n ": 


'o ' ' 

^ : 


ma : 

. 1 .... 1 .... 1 .... 1 . 


:m 

. 1 .... 1 .... 1 .... 1 . 


. 1 .... 1 .... 1 .... 1 . 


:'Vl7.b' 

\ <m> : 

. 1 .... 1 .... 1 .... 1 . 


:'Vl8.5 : 

<% 

. 1 .... 1 .... 1 .... 1 . 





00" 

-05°35'10" 



18 h 38 m 41.5 s 40.5 s 




R.A. (J2000) 

Fig. 6. — SMA 13 CO (2-1) velocity channel maps. For channels of 96-112.5 kms , the 
maps made with the combined SMA and IRAM 30m data are shown on the bottom. The 
contour levels start from and increase in steps of 0.3 Jybeam -1 . 
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Fig. 7. — (a) A close-up of the SMA CO (2-1) velocity integrated map. Blue and red 
contours show emission integrated from 69 to 100.5 kms" 1 and from 106.5 to 139.5 kms -1 , 
respectively. The contour levels start from 20% and increase in steps of 10% of the peak 
emission, which is 59.1 and 71.3 Jybeam -1 kms -1 for the blue and red lobes, respectively, 
(b) The same as (a), but for 13 CO (2-1), showing emission integrated from 85 to 102 kms -1 
and from 105 to 120 kms -1 . The peak emission is 37.7 and 26.2 Jybeam -1 kms -1 for the 
blue and red lobes, respectively. 
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Fig. 8. — (a): The color image shows the first moment map of the CH3OH (15^11-163^3) E 
emission. The gray contours show the zeroth moment map starting from 1.0 Jykms -1 and 
increasing in steps of 0.8 Jykms -1 . A dotted line marks the P-V cut at PA=140°. Two 
arrows are the same as those in Figure [HI indicating the outflow orientation, (b)-(c) Same 
as (a), but for CH3CN (I27-II7) and vibrational CH3CN (12i-lli). The contour levels in 
panel (b) start from 0.6 Jykms -1 and increase in steps of 0.8 Jykms -1 , and in panel (c) 
start from 1.3 Jykms -1 and increase in steps of 0.7 Jykms -1 . The color wedge on the right 
of panel (c) indicates the velocity scale of the first moment maps, (d)-(f) P-V diagrams 
shown in both gray scale and contours. The contour levels all start from 0.21 Jybeam -1 , and 
increase in steps of 0.35, 0.21, and 0.28 Jybeam" 1 for CH 3 OH (15 4 ,ii-16 3 ,i3)£ (d), CH 3 CN 
(I27-II7) (e), and vibrational CH 3 CN (12i-lli) (f), respectively. The horizontal line marks 
the peak position of MM1, which is taken to be the zero offset. The vertical line denotes the 
systemic velocity. 
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Fig. 9. — (a) Solid and dashed contours show the emission integrated from 96.6 to 102.0 
kms" 1 and the absorption integrated from 103.2 to 106.8 kms -1 , respectively, in CN (N=2- 
1, J— 3/2-1/2, F— 5/2-3/2). The starting and spacing contour levels are 54 mJykms -1 for 
the emission and —63 mJykms -1 for the absorption, (b) The spectrum of the same CN 
transition, extracted from MM1. A vertical line marks the systemic velocity. The velocity 
channels integrated to show the emission and absorption in (a) are shaded with slashes and 
dots, respectively. 
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Fig. 10. — The 13 CO (2-1) P-V diagram in gray scale, overlaid with those in C 18 (2-1) 
(a) and CH3CN (I27-II7) (b) in contours. The C 18 starting and spacing contour levels 
are 0.21 Jybeam -1 . In each panel, the red solid curve shows the model of a free-falling and 
Keprelian-like rotating envelope, and the blue dashed curve represents the same model but 
with scaled down parameters (see Section fl~2l for more details of the models). 



Table 1. Summary of SMA Observations 







Freq. Cove. (GHz) 












Obs. Date 


Array Conf. 


LSB 


USB 


Spec. Res. (kHz) 


Bandpass Cal. 


Gain Cal. 


Flux Cal. 


T"225GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


2007 May 17 


VEX (6) 


219.3 - 221.3 


229.3 - 231.3 


406.25 


3C273 


J1743-038, J1911-201 


Vesta 


0.1 


2007 Jul. 08 


Comp. (8) 


219.3 - 221.3 


229.3 - 231.3 


406.25 


3C273 


J1743-038, J1911-201 


Vesta 


0.07 


2008 May 21 


Comp. (8) 


215.4 - 217.4 


225.4 - 227.4 


812.5 


3C273 


J1743-038, J1911-201 


Callisto 


0.15 


2008 Jul. 19 


EXT (6) 


215.4 - 217.4 


225.4 - 227.4 


406.25 


3C273, 3C454.3 


J1743-038, J1911-201 


Uranus 


0.07 



Note. — Col.(l): Observing dates. Col. (2): Array configurations, with the Very Extended, Compact, and Extended configurations abbreviated to VEX, 
Comp., and EXT, respectively; the figure in parentheses denotes the number of antennas used for each observation. Col. (3)— (4): Approximate ranges of 
rest frequencies covered in the lower sideband (LSB) and upper sideband (USB). Col. (5): Spectral resolutions. Col.(6)-(8): Bandpass, time dependent gain, 
and absolute flux calibrators. Col. (9): The averaged 225 GHz opacity. 



